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Abstract 

We study the magnetic instability of gapless superconductors. The instability arises due to the 
infrared divergence associated with the gapless modes of the superconductor. When the Fermi- 
surface mismatch between pairing fermions is close to the gap, the gapless modes have a quadratic 
energy dispersion relation at low energy and open a secondary gap at the Fermi surface, which 
is only power-suppressed by the coupling. On the other hand, for a large mismatch, we find the 
gapless superconductor does not open a secondary gap, but instead makes transition to a new 
phase by forming the condensate of supercurrents. We calculate the condensate of supercurrents 
by minimizing the effective potential. In the new phase, the Meissner mass is positive for the 
magnetic fields orthogonal to the direction of the condensate but zero in the parallel direction. 
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How matter behaves at extreme density is ultimately related to the fundamental questions 
on basic building blocks of matter. Quantum chromodynamics (QCD), which is believed 
to be the theory of strong interaction of elementary particles, predicts matter at extreme 
density is color superconducting quark matter, where the quarks form Cooper-pairs 1]. The 
search for quark matter in heavy ion collision or in compact stars is currently under intense 
investigation. 

Quark matter is expected to exhibit a rich phase structure, having various pairing patterns 
for quarks, as temperature or density changes. At high density, where the SU(3) flavor 
symmetry among u, d, s quarks is good, the QCD interaction pairs quarks to form color- 
flavor- locked (CFL) matter, respecting the flavor symmetry [2j. The flavor-asymmetric 
electroweak interaction and quark mass strain the CFL matter, which then becomes unstable 
under a large stress. 

BCS pairing breaks in general when the stress is bigger than the pairing energy or 8fx > 
Ao/\/2, where the stress 2 8 fx is the chemical potential difference of pairing quarks and Ao 
is the BCS gap Q. Recently, however, it has been shown that Cooper-pairing may be stable 
even if the stress is bigger than the gap, known as Sarma phase |4], if one enforces the electric 
neutrality in quark matter or uses a momentum-dependent kernel for pairing The 
salient feature of such asymmetric quark matter is that quarks can be excited at arbitrarily 
low energy, though it is superconducting. The gapless superconductors, if realized, will be a 
new kind of Fermi liquid, which has both properties of the usual BCS superconductor and 
the Landau-Fermi liquid. It was soon found, however, that the gapless superconductors have 
negative Meissner mass squared |7| or a negative supercurrent density .8], showing magnetic 
instability of the system. 

In this letter we calculate the effective potential for the gapless superconductors and 
study the magnetic instability. We find that the magnetic instability arises because of the 
infrared divergences associated with the gapless modes. When the stress is greater than the 
gap, 8 fx ^ A, the infrared divergence is logarithmic and the system develops a condensate 
of supercurrents. On the other hand, when the stress is very close to the gap, 8 fx ~ A, the 
divergence is much more severe and it inevitably leads to opening a secondary gap at the 
Fermi surface. 

To study the secondary gap in detail, we derive the low energy effective Lagrangian for 
the gapless modes by integrating out the gapped modes, which turns out to have an attrac- 



2 



tive four-Fermi interaction. Unlike the ordinary Fermi liquid, the four-Fermi interaction in 
gapless superconductors scales like s~ l l 2 as one scales down to the Fermi surface (s — > + ), 
if incoming fermions have equal and opposite momenta [10], and induces a gap which is 
only power-suppressed in couplings. Both renormalization group (RG) analysis and gap 
equation analysis show that the gapless excitations develop a gap, which naturally stabilizes 
the system. 

Let us consider a minimal model for gapless superconductivity, which has two different 
flavors that pair, denoted as up and down, and free electrons for electric neutrality, shown in 
Table |U The pairing force is SU(2) C color interaction under which the up and down particles 
are fundamental. The system is not a color superconductor, since the Cooper pair is a color- 





chemical potential 


electric charge 


SU(2) C 


V>i (up) 


fii = p> — Sn 


qi = q + 5q 


2 


ip2 (down) 


fJ-2 = + 5[1 


q 2 = q- 5q 


2 


tp e (electron) 


He = 25fl 


Qe = -1 


1 



TABLE I: The electron chemical potential is // e = — and q e = —25q, since the system is in 
equilibrium under the weak interaction, ip2 <-> ipi + e ~ + ^e- 

singlet. However, it is an electric superconductor and exhibits all the essential features of 
gapless superconductivity, including the magnetic instability. It is straightforward to extend 
to the gapless 2-flavor color-superconductor (g2SC) or to the gapless CFL superconductor 
(gCFL). Neglecting the anti-particles, the system is described by a Lagrangian density 

£ = Z>J [id t - E(p) + m] + ! e^lV^'Vi'Vy + CL , (i) 

1=1,2 

where E(p) is a flavor-independent function of momentum and the color and spin indices 
are suppressed. Gluons are integrated out to generate an attractive four-Fermi interaction, 
antisymmetric in color and flavor, and also other irrelevant interactions, denoted as £' nt , 
assuming the gluon exchange interaction is most attractive for quarks antisymmetric in 
color and flavor. 

Neglecting the irrelevant interactions, one can calculate the free energy in the mean field 
approximation, 

A 2 f d 3 v 
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Ve 2 + A 2 + Sfi + Ve 2 + A 2 - Sfj, 



12tt 2 



(2) 



where e(p) = E(p) — fx and the last term is the electron free energy. When 8 fx > A, the 
gap equation has an unstable solution, A = a/ A (28 fx — A ), where A is the gap when 
8 fx = 0. This unstable Sarma phase can be stabilized if charge neutrality is enforced. The 
electric charge neutrality condition for the Sarma phase, ((Jem) — 0, is satisfied for fx ^> 8 fx 
if 8fx = + A 2 , where 2 8fXo is the chemical potential difference of the free neutral 

system Compared with the free neutral system, the free energy becomes for fx ^> 8 fx 



fx 2 

tt s (A, 8fx) - Q {ree ~ — - 8fx 2 

7T Z 



(3) 



The neutrality enforcement stabilizes the Sarma phase if A /4 < 8fx Q < A /2 or the gap 
< A < A /v2- At the minimal 8fi , the ratio 8fx/A is smallest and close to 7/4 — l/\/2. 

We now calculate the Coleman- Weinberg effective potential [111] for photon fields in the 
Sarma phase, which will be same as the free energy in Eq. (j2J) except that p — > p + qi eA 
and p — > p + qi eA As we integrate out the high frequency modes, the effective potential 
will get non-local terms, which can be expanded in powers of momentum. When all the 
modes with u > 8 fx are integrated out, only the gapless modes survive and the effective 
potential becomes 

2" 



V(A) = ~ 4 Jj2^ ln Pa + { veie2 + A 2 - 8 fx 



+ C (A) + ^^A--^^A 2 + - ■ ■ (4) 



where 8fx = 8 fx — 8qeA and = E(p + q,ieA) — fx — qeA . In the effective potential we 
have introduced an ultra-violet cutoff A (< 8 fx) and the counter-terms, Cj's and B^s. The 
ellipsis denotes the higher order terms in photon fields and their derivatives. 

As we will see, the gapless superconductor suffers instability due to the infrared diver- 
gences associated with the gapless quasiparticle modes. Depending on how close the stress 
is to the gap, the structure of divergences and their physical consequences differ drasti- 
cally. When the stress is not too close to the gap (8 fx > A), as in g2SC, the quasi-particles 
have approximately linear dispersion relation near the Fermi surface. In this case there is 
a logarithmic infrared divergence and the system is stabilized by spontaneously generating 
Goldstone currents. However, when the stress is close to the gap (8 fx m A), as in gCFL, the 
quasiparticles near the Fermi surface have approximately quadratic dispersion relation and 
there is genuine instability in the system, which leads to opening a secondary gap at the 
Fermi surface. 
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We consider first the case where the stress is not too close to the gap. If we integrate out 
further the quasi-particle modes till the remaining has an approximately linear dispersion 
relation, the effective potential becomes 



V(A) 



47T 



(eVi • A f 



In 



A 2 



{eV 1 ■ A)' 



+ 1 



+ (1 



(5) 



+ ^o(A) + ^ 2 -^V + 



2 2 

where V\ and v% are the density of states at the Fermi surfaces at p\ and p2, respectively 
At the inner and outer Fermi surfaces, e(pi) = — a/5/U 2 — A 2 and efa) = \J — A 2 . 
The angular integrations are over the inner and outer Fermi velocities, v% = dE/dpi and 
V2 = 0E/dp2- We have also defined rj = 5p,/ ^ bp? — A 2 , Vf = (rjq + 5q,rjqvi), and V£ = 
(rjq - 5q,r]qv 2 ). 

We see that the one-loop effective potential due to the gapless modes is negative for both 
Aq and A, and its second derivative is negative infinity at the origin. Since the Meissner 
mass term in the potential has the positive sign, while the Debye mass has a wrong sign, 
the effective potential has a minimum away from the origin for A. We fix the counter-terms 
by imposing the renormalization conditions at a scale M, 



d 2 V 



PlA 2 

un -0 A =M 

to get the effective potential 

1 



-m 



D i 



1 d 2 V 
3 ij OAdA, 



m 



M ? 



(6) 
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M 2 
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^m 2 D Al --{vx + v 2 ) {m + 5qf e 2 A 2 

7] 
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(7) 



The physics lies in the renormalization conditions that we imposed in Eq. (JHJ). If we take 
M = 0, we will recover the negative mass squared for the Meissner mass, similar to the result 
obtained in 13]. However, this result is sensitive to the ultraviolet cutoff. Here instead we 
take a nonzero M such that the Meissner mass due to modes with uj > M is nonnegative 
and then we study how the system flows as we change M, keeping the ultraviolet physics in 
the renormalization conditions. 

The vector fields in Eq. (J7J) are not physical, because the effective potential we derived 
is not manifestly gauge- invariant. For the physical degrees of freedom, we go to a unitary 
gauge, where the vector fields are replaced by the gauge-covariant combination of gauge 



fields and Goldstone fields, A — > A — Vy>, A — > A + <9 t y2. By minimizing the effective 
potential, Eq. (JJJ), we find that the system develops a condensate, breaking the rotational 
invariance, 



A — Vip) 2 )=M 2 exp 



2rj [y\v\ + V2V2) e 2 q 2 

Since the condensate is invariant under the renormalization group flow, we calculate it 
at M = 5fi. The Meissner mass due to the fast modes (00 > 5 fx) is found to be m 2 M ~ 
(8/3) q 2 e 2 v 2 , where and ti* are the density of states and the quasiparticle velocity at 
the pairing momentum, p*, respectively. We find 



(a — V(p ) / — ^a* 2 ex P 



77 [y x v\ + v 2 v\) 



(9) 



- .=011 

A=A C 



which shows the system has a spontaneously generated Goldstone current even in the absence 
of external gauge fields, ^— V<^ = A c 7^ jl4j . The ground state, however, does not carry 
any net current, since = ~~ dV(A)/dA 

The gapless superconductor is a directionally perfect diamagnet because the effective 
potential has flat directions along the angular rotation of A c . Under an external field A, the 
current is given as 



Ji{x) 



D 2 V(A) 



(A j - d jV ) . (10) 

A=A C 



dAtdAj 

Taking the curl of the current, we get V 2 B = c 2 M (b — B ■ A c A^j, where where A c is the 
unit vector along the condensate and c 2 M = ^ {^xv\ + V2V2) e 2 q 2 . The Meissner mass is cm 
for the magnetic fields orthogonal to the condensate and zero for the parallel fields. 

When 5fi ~ A, pi « p 2 ~ and the gapless modes have a quadratic dispersion relation 
u(p) ~ (v* ■ /) 2 /(25/i) for < to < 5fi, where the residual momentum I = p — p*. The 
effective potential due to the quadratic gapless modes is found to be, with K = (1, if*), 

^ (A)= _ ^(5/4)r ( i/4) ^^ W 2 l*v..^ + W 2 . (ID 
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The effective potential shows a genuine instability, not stabilized by the condensate of Gold- 
stone currents, since it is unbounded from below. However, as we will see later, by the 
Kohn-Luttinger theorem Q], the quadratic gapless modes pair among themselves to open 
a gap at p — p* and the system becomes stable. 

The Sarma phase has a gapless excitation, denoted as whose energy spectrum is 
given as to = ± (d/j, — a/ e(p) 2 + A 2 j , and a gapped excitation, denoted as \1/h, with energy 



spectrum uo = ± (S[i + a/ e(jo) 2 + A 2 j . The energy of the gapless mode, measured from the 
Fermi surface is quadratic in the residual momentum except very near the Fermi surface. 

Subtracting out the Fermi momentum as in the high density effective theory 17], we may 
write the effective Lagrangian for gapless modes as, when 5fi > u > u;ir = (<5/i 2 — A 2 ) / {25 ji), 



eff 



id t + 



$?(v*,x) + -\l> t \l/ t v|>\l/ H 



(12) 



where the summation is over all the patches that cover the Fermi surface and the ellipsis 
denotes the higher order interactions. 

Because of the quadratic dispersion relation, as we scale down to the Fermi surface uo — > 
s oo (0 < s < 1), the momentum parallel to the Fermi velocity v* scales as u* • I — > s 1 / 2 v* ■ I, 
while the perpendicular momentum l± — I — v* v* • I does not scale. Since the action for 
the kinetic term has to be scale-invariant, the gapless mode scales in the momentum space 
as *&(t, I) — > s™ 1//4 ^(t, I ) . The immediate consequence of such unusual scaling is that 
the four-Fermi interaction becomes a relevant operator when the incoming fermions have 
equal and opposite Fermi momenta. Under the scale transformation the four- Fermi coupling 
transforms as k — » s -1 / 2 k and it will hit an infrared singularity unless a gap opens in the 
infrared region. 

Among the irrelevant interactions in C' int , which we have neglected so far, there is a 
repulsive four- Fermi interaction, C[ nt 3 Gr a ViV'i'02'02> which is symmetric in color. This 
interaction describes transition between the gapless mode and the gapped mode and induces 
an attractive four-Fermi interaction for the gapless modes. In our minimal model, the 
attractive channel turns out to be a spin-1 and color-singlet channel and the condensate 



takes (\I / zcr2cr 4 A 2, I / ) ~ A s <5* 3 , similar to the polar phase Q], where cr's are the Pauli matrices 
in the spin space and A 2 is the antisymmetric Pauli matrix in the color space. By integrating 
out the gapped modes, we have k = (7r/4) G 2 s v+Jv* at one-loop (See Fig. 1). 
The Cooper-pair gap equation is given as 

d 4 / A s 



A., 



(2k) 



P 



1 2 



Upon integration, we get 



6.85 



25/i 



4.2 



(13) 



(14) 
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FIG. 1: The four- Fermi interaction of the gapless modes induced by the gapped modes VPjj. 

where g' 1 = In (2/Z/A) . For a weak coupling, g < 1, the secondary gap is well separated 
from the ultraviolet scale. For a strong coupling, however, one needs to go beyond the mean 
field approximation to find the correct secondary gap. In any case, the gapless supercon- 
ductors are stabilized by opening a secondary gap at the Fermi surface if the secondary gap 
is bigger than the infrared cutoff of the effective theory, u>ir — (5fi 2 — A 2 )/ (25 n). 

We have studied a minimal model for gapless superconductivity and found that the gapless 
superconductors necessarily open a secondary gap at the Fermi surface if the characteristic 
scale of the gapless superconductors cufr is small enough compared to the secondary gap. 
However, if uir > A s , the gapless superconductors do not open a gap but instead develop 
a condensate of supercurrents. 

Though our analysis is generic and applies to any gapless superconductors or superfluids, 
the specific form of the secondary gap depends on the details of the gapless superconductors. 
For instance, for g2SC the gapless modes are degenerate in color and spin. Thus, if the 
secondary gap forms, it will be a color-antitriplet and spin-1 gap, but not exponentially 
suppressed. For gCFL the gap has to open in the color-sextet channel, since the gapless 
modes are not deg enerate in color. The possible candidate is a color-sextet condensate, 



To conclude, we have shown that the gapless superconductors are unstable due to the 
infrared divergence associated with the gapless modes and make phase transition either 
to superconductors with a supercurrent condensate, or to gapful superconductors. The 
Meissner mass is nonnegative in both cases. In the former case, the rotational symmetry is 
spontaneously broken and the Meissner mass is directional. In the latter case the secondary 
gap is not exponentially suppressed but only power-suppressed in couplings, which may have 
significance in neutron stars or in atomic superfluids. 



(^lCjq^r) 19]. 
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